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Abstract: It has been proposed that during the formation of the first generation stars 
there might be a "dark star" phase in which the power of the star comes from dark mat- 
ter annihilation. The adiabatic contraction process to form the dark star would result 
in a highly concentrated density profile of the host halo at the same time, which may 
give enhanced indirect detection signals of dark matter. In this work we investigate the 
extragalactic 7-ray background from dark matter annihilation with such a dark star for- 
mation scenario, and employ the isotropic 7-ray data from Fermi-LAT to constrain the 
model parameters of dark matter. The results suffer from large uncertainties of both the 
formation rate of the first generation stars and the subsequent evolution effects of the host 
halos of the dark stars. We find, in the most optimistic case for 7-ray production via dark 
matter annihilation, the expected extragalactic 7-ray flux will be enhanced by 1 — 2 orders 
of magnitude. In such a case, the annihilation cross section of the supersymmetric dark 
matter can be constrained to the thermal production level, and the leptonic dark matter 
model which is proposed to explain the positron/electron excesses can be well excluded. 
Conversely, if the positron/electron excesses are of a dark matter annihilation origin, then 
the early Universe environment is such that no dark star can form. 
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1. Introduction 

It has been proposed that there might be a new phase of the stellar evolution, dark (matter 
powered) star (DS), during the formation of the first generation stars (or Pop III stars) in 
the early time of the Universe (e.g., [|l], |2|, |^, ^, ^, 0]). Collapse of the gas would attract 
dark matter (DM) to collapse together and form a high density core. Self-annihilation of 
DM particles would inject energy into the gas inside the core, which heats the gas core and 
prevents further collapse of the gas to initiate nuclear fusion. DSs typically have distinct 
features compared with normal stars e.g., they are massive (500 — 1000 M0), large 
(1 — 10 AU), and have low temperature (Tguj-f < 10000 K). They might be detected with 



large telescopes j^, |10| 



The impact of DS formation is very important in many aspects of astrophysics and 
cosmology. The reionization history of the Universe would be altered due to different 



contribution to the ionization photons from DSs [11|. The DS remnants might provide 
the seeds for super-massive black holes [||]. Moreover, the formation of DSs would result 
in a cuspier density profile of DM outside the stellar core^, which is expected to give an 
enhanced annihilation luminosity and can be reflected in the indirect detection signals such 
as 7-rays and neutrinos |]ll|, |l^] . 

In this work we aim to discuss the contribution of DM annihilation to the extragalactic 
7-ray background (EGRB), within the framework of DS formation models. Schleicher et al. 
have studied the extragalactic background radiation for the DS scenario and derived the 



constraints on DM particle parameters |11]. There are several improvements or updates 



in the current study: 1) A more realistic prescription on the formation of Pop III stars 



^Strictly speaking, it is the adiabatic contraction which results in a cuspier DM density profile, and due 
to the same reason a DS is formed. In the following when we say "DS enhancement" it actually means the 
enhancement induced by the contraction process to form DS. 
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(which can give birth to the DSs) is adopted; 2) The DM annihilation models are discussed 
extensively, especially for the leptonic DM models which may be responsible for the recent 
reported electron/positron excesses m, HI, P, H, 0; 3) The EGRB data is updated with 



the new result released by the Fermi team |18], which gives stronger constraints than the 
previous EGRET data 

Throughout the paper a flat ACDM cosmological model is adopted, and the cosmo- 
logical parameters are adopted as the results of the combined analysis of the WMAP 5 
year data and other cosmological measurements: 0,m = 0.274, i}\ = 1 — Oa/, = 0.046, 
h = 0.705, as = 0.812, Ug = 0.96 [|^]. Adopting a different set of cosmological param- 
eters does not alter the final results significantly. For example we have checked that the 
formation rate of first stars would change by less than several percent if we adopt the 
cosmological parameters derived with the WMAP 7 year data together with the baryon 



acoustic oscillation and Hq data |21]. 



This paper is organized as follows. In Sect. 2, we present the formation rate of first 
generation stars, which is tightly connected with the formation of DSs. In Sect. 3, we 
calculate the annihilation luminosity and 7-ray flux of both the scenarios with and without 
DS formation, and derive the constraints on DM model parameters. The conclusion and 
discussion are given in Sect. 4. 

2. Formation rate of the first generation stars 

We first introduce the mass function of the DM halos. The comoving number density 
distribution of DM halos can be expressed as 



where v = 5c{z)/cr{M), 6c{z) = 1. 68/1^(2;) is the critical over-density in spherical collapse 
model, D(z) is the linear growth factor |22[ |. A, a and p are constants. For (A, a, p) = 
(0.5, 1, 0) it is the Press-Schechter (PS) formula @, and for {A, a, p) = (0.322, 0.707, 0.3) 
it is Sheth-Tormen (ST) formula |2j]. In this work we adopt the ST mass function. a^{M) 
is the average variance of density field 



a\M) = ^ J W\kRM)P5{k)k^dk, (2.2) 



with the top-hat window function W{x) = 3(sin3; — xcosx)/x^. Ru = (3M/47r/9m)^^^ is 
a radius within which a mass M is contained with a uniform matter density field. The 
matter power spectrum Ps{k) is expressed as 

Ps{k)=As{k-Mpcr^T\k), (2.3) 

where Ag is normalized using erg, transfer function T{k) is obtained from a fit of CDM 
model [25 1, 

Tig) = ^"^^+3^^^^^^^ [1 + 3.89, + (16.1,)^ + (5.46,)^ + (6.71g)^] (2.4) 
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with q = k/hV and F = Q.mh eyi.\i[—VL}y{l + ^/2h/VLm)\■ 

In the early Universe, the first stars would form in halos that are massive enough 
for gas to cool and condense. The number of such halos increases rapidly, and the halo 
destruction by mergers can be neglected. We therefore assume that the Star Formation 
Rate (SFR) of the first generation stars is proportional to the redshift derivative of the 



halo mass function within a certain mass range [26|, with a time delay Td due to the 



cooling and collapse of the halo |27]. At a given redshift z, the SFR can be written as 

-MmaxC^") 



SFRpopiii(^) 



Afmin(^") 



[1 - Qh+(^")] ^[z - z'{Td, z")]dM'd/', 



dMdz 

(2.5) 

where Qh+ is the volume-filling factor of H^" regions, z'{Td,z") is the redshift with time 
Td delay of z" . Here we further assume that only one first generation star could form in 
each potential halo'^. The 6 function in Eq. ( |2.5| ) means that the halos which satisfy the 
formation condition of the first stars will eventually contribute to the SFR Td later. For the 
time delay, we adopt Td = TcooI + ''"ff j in which the cooling time scale for H2 is TcooI = 

2.38 X 

10^^ (10^) '^ '^^^ i^y^'^^ s and the free fah time scale is rg = 2.77 x 10^^ (^)"^^^ s 
p7|]. This delay time scale Td is mass-dependent, which makes the calculation of Eq. (p 



more complicated. 

As for the lower mass limit, it is usually believed that the first generation stars could 
only form in halos with virial temperature above 10^ K if the coolant is H2. However, 
in the presence of Lyman- Werner (LW) photons which are emitted by stars which form 
previously, this lower limit increases since only H2 in massive enough halos can survive from 
the photo-dissociation. Following Ref. [^^, the lower limit of the halo mass Mmin for the 
first star formation would be the maximum of Mhj-cooI = 6.44 x W^MqJ2i'^'^ (^sf) 
and Mvir(10^i^, z), with J21 the specific intensity of LW radiation in unit of 10~^^ erg 
s~^cm~^Hz~^sr~^. 

The upper limit can be adopted as the virial mass corresponding to a virial temperature 
of 10^ K. It is generally believed that within halos with Tvir > lO'* K, the cooling by atomic 
Hydrogen excitation becomes more efficient, so that the gas would fragment into multiple 
parts and the massive and isolated first generation stars cannot form. 

The intensity of LW radiation from the first generation stars is 

npopm{z'){l + z'feiMj^dz', (2.6) 

in which eLw is the emissivity of LW photons, ^ is the comoving distance per redshift, 
Zs is the upper limit above which the LW photons would be redshifted out of this band, 
i.e. Zs = 13.6/11.2(1 -|- z) — 1. npopiii(z) is the comoving number density of active first 
generation stars, which can be obtained through integrating SFR with respect to redshift 
within the life time of such stars 

"-Popiii(2) = / SFRpopiii(z')d2'- (2-7) 

J Z + Az{TpopUl) 

^There has been no discussion on whether DSs can form if fragmentation is important during first star 
formation (e.g., |28|). Our discussion does not apply to those scenarios. 
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We adopt a typical mass of 200 Mq of the first stars in this work. According to Ref. |]29| , 
the life time of such a star is ~ 2.24 million years, and eLw — 1-15 x 10^^ erg s~^Hz~^. The 
choice of this mass is consistent with the numerical simulations of first stars. If the mass 
parameter alters by a factor of several, the main infiuence would be the LW feedback. As 
discussed below, we will investigate the LW feedback efficiency within a wide range. This 
effectively includes the variance of the first star mass to some extent. We also note that if 
the initial mass function (IMF) of the first stars follows a Salpeter form between 100 Mq 
and 500 Mq [3C], the LW emission rate per unit star mass is close to that of a 200 Mq 
star. 

Note that DS formation may also affect the LW feedback. First, DSs will also produce 
LW photons. However, for typical DS parameters, L ~ 10^ - 10^ Mq, T ~ 5000 - 10000 K 
Q, we find that the contribution to LW emissivity (11.2 — 13.6 eV) is negligible. Second, 
after the DS stage, the object will enter the main sequence as a traditional star powered 
by nuclear fusion, with different properties from the ordinary PopIII star. There is no 
clear conclusion about the properties and fate of this particular main sequence star (e.g., 
||8|, 1^ |3l|, 0]). In any case, since we have employed a large range of uncertainty of LW 
feedback (see below), the impact of DS formation on the first star formation rate would 
not introduce much larger uncertainties, and we have neglected this effect in the following 
discussion. 

Besides the first generation stars, there are also galaxies that can contribute to the LW 
dissociation and ionization process. SFR of metal-enriched and metal-free galaxies are: 

SFR-1 = pen(^) X -^/. / JirriM, z)dM (2.8) 



and 



SFR%, = [1 - pen(^)] X -^A / TJ71-iM, z)dM (2.9) 

respectively, where = 2 x 10^'^ is the fraction of gas that converts to stars, and pcn is 
the enrichment probability of gas by metals. The enrichment probability could be divided 
into two parts: the self-enrichment through which a Pop III star could pollute the host 
halo, and the probability that a halo is polluted by metals from neighbors. We follow Ref. 
1 33] to calculate the self-enrichment probability. The pollution probability from neighbors 
is assumed to be just the metal fraction of the Universe. According to Ref. almost 
one half of the mass of a Pop III star would be released through pair instability supernova 
(PISN), and the metal yield from galaxies is about 0.04 As for the LW photons 



emitted by galaxies, we take a time-evolved emissivity from Ref. [36|, and add this part to 
Eq. (P). 

We also trace the reionization evolution history, since it could reduce the SFR of first 
generation stars by a factor 1 — Q-^+ . The evolution of (5h+ is described as [ pT| ] 

dQH+ _ Qu+C(,z)n^+aA driph/d^ 
dz H{z){l + z) nn ' ^ ' 

where C{z) = 27.466 exp(-0.114z + 0.001328z^) is the clumping factor is the 



number density of ionized hydrogen, oa is the recombination coefficient ||3^, H{z) is the 
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Hubble parameter, and nn is the mean neutral hydrogen density. In our work, the produc- 
tion of ionization photons dnp\i/dz is expressed as 



^(.) = (! + .) 



npopiii(^;)eph + / SFRgai(2;')egai(2; - z')dz' 

J oo 



dt 



The first term refers to the contribution from first generation stars. For a 200 Mq metal- 
free star the emissivity of ionization photons eph is 3.1 x 10^" s~^ [^]. The second term 
is the contribution from the galaxies, including metal-enriched and metal-free galaxies. In 
the second term, egai is a time-evolved emissivity of ionization photons, which is adopted 
from Ref. ||. 

There are still large uncertainties about the LW feedback mechanism due to the lack of 
understanding of the self-shielding pg| ]. Furthermore, the molecular hydrogen in the relic 
HII regions may also absorb some LW photons, reducing the effect of photo-dissociation. 
To consider such a reduction on the LW intensity, we employ a free parameter /lw to 
represent the efficiency of LW feedback. 



Combining Eqs. (|2.5[)-( 2.11 ) together, we can get a coupled equation set of the SFR 
at a given redshift. This equation set depends on the evolution history of SFR, so it should 
be integrated step by step from an initial time. At a high redshift, e.g. around 60, Jlw 
is negligible, SFR could be obtained directly by setting Mmin = Mvir(10^K, z). With this 
initial value, the SFR at first step after the initial redshift could be calculated by the set 
of equations discussed above. Performing the same calculation step by step, we can finally 
obtain the entire evolution history of the SFR. 

The SFRs of Pop III stars for different LW efficiencies /lw = 0.0, 10"^, IQ-^, and 1.0 
are shown in Fig. |l|. It is shown that for the very effective LW feedback the star formation 
ends very early, at z ~ 25. For the most optimistic case without LW feedback the star 
formation can last till z ~ 6. These two extreme cases may represent the lower and upper 
limits of the SFR of the first generation stars. 



3. EGRB from DM annihilation 

3.1 Enhancement factor of DM dumpiness 

In this section we describe the enhancement factor of DM annihilation due to dumpiness, 
defined as the ratio between DM annihilation rate with structures and that of smoothly 
distributed case. We first consider the standard case without DS formation. The details 
of the calculation of the enhancement factor of the extragalactic DM annihilation can be 
referred to Refs. [^, Here we quote the basic formulae. The density profile of all the 
halos is assumed to be NFW profile [^3| 

Pir) = I , — ^-yi (3-1) 

{rlrs){\ + rlrsY 

with two scale parameters and ps . The scale radius is determined by the concentration 
parameter of the halo 
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Figure 1: The evolution of SFR of first generation stars. In this plot the SFR is the result of Eq. 
(2.5) multiplied by Affgdz/di to convert to the usual definition. The solid curve is the fit to the 
observational SFR for the second and third generation stars . 



where r„ is the virial radius. The scale density ps is then determined by normalizing the 
mass of the halo to M. Note that there is no consensus about the density profile at the 



innermost part of the halo ||4^, 45, 46, 47, For different density profiles of the halos, the 
annihilation signals would differ from each other very much. However, as shown in [||], the 
density profiles of the DM halo after adiabatic contraction are very similar to each other 
even for very distinct initial density profiles. Since in this work we mainly focus on the 
effects of the adiabatic contraction (DS formation), we do not discuss in detail the effect 
of different halo density profiles. 

The concentration parameter as a function of halo mass M can be extracted from 
N-body simulations. For the low mass halos which are beyond the resolution of simu- 
lations, the concentration is obtained by extrapolation. In this paper we consider two 



concentration models. One is the semi-analytical model developed in Ref. [49| with the 
update of WMAP5 cosmological parameters, which is labeled as "BOl". The other model 
is to extrapolate the power-law fitting results from the numerical simulations based on the 
WMAP5 cosmological parameters ||5^ directly to the low mass range, which is labeled as 
"power-law". For both models we employ the redshift evolution Cv{z) = Cy{z = 0)/(l -|- z) 



|49]. The results of BOl and power-law models do not differ very much from each other 
at high mass scales (M > 10^ Mq). However, when extrapolating to lower mass scales, 
power-law concentration would be larger than BOl model, and hence gives larger annihi- 
lation signals. The halos which can give birth to DSs are massive enough that the DS 
enhancement between these two concentration models will be very similar. There are some 
updated concentration models which show different mass-dependence and redshift evolu- 
tion (e.g., ^]), however, they will also suffer from the uncertainties when extrapolating 
to low mass halos. 

The total annihilation luminosity for a population of DM halos with comoving number 
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density distribution dn(z)/dM is given as 

L,ot{z) = I dM^(l + zfL{M, z), (3.3) 

where (1 + z)^ is to convert the comoving halo mass function to the physical one, and 
L{M,z) = An J^" p^{r)r'^dr is the annihilation luminosity of a single halo. For the mass 
function dn(z)/dM please refer to the Appendix. The enhancement factor is then 

where is the average DM density of the Universe today. 
3.2 DM annihilation luminosity vi^ith DS formation 

If there is DS formation following the first generation stars, the annihilation luminosity 
will be enhanced. The annihilation luminosity of the halos which ever contained DSs can 
be calculated as 



Lt,{z) = {l + zf r dz' 

X / / UMM,z"')[l-Q^,{z"')]-^^{M,z"') 
Jz' Jm^^Jz'") dMdz' 



Afn.in(^"') 

X 6[z' - z"{Td, z"')]dMdz"', (3.5) 

where fds is the enhancement factor of the annihilation luminosity of a single halo after 
DS formation compared with the original density profile, and r^s is the age of the halo 
which hosts a DS. Generally speaking the density profile would steepen as r~^'^, regardless 
of the initial density profile of the DM halo. We note, however, that the detailed density 
profile would depend on the parameters of the initial DM halo and the collapsing gas core. 
According to the density profile of halos after DS formation [p^], we estimate fds ~ 10^ 
compared with initial NFW profile. Such a value is consistent with that estimated in 



If there is only an annihilation effect to consume the DM, t^s is estimated to be much 
longer than the age of the Universe, and the lower limit of the redshift integral {dz') can be 
taken as oo. On the other hand, in order to qualitatively take into account the evolution 
effect of the halos after the DS formation, we also take another conservative approach 
by considering that the halo with DS formation can only exist for a finite time due to 
possible mergers with other halos^. Following Ref. ||3^, for a halo with mass M ~ 10^ 
M0 at redshift Zi, we define a surving time Tds (corresponding to a final redshift Zf) after 
which the halo reaches a mass 2M due to merger or accretion. There will be a probability 



^According to the numerical simulations of Ref. jH^, the result of a merger between two halos with 
different density profiles could attain a density profile with an intermediate cusp slope. However, in Ref. 
JssI ] it was found that the final density profile would be essentially close to the cuspier one. Here we assume 
the initial density cusp r^^'' of the dark-star-enhanced halos gradually shallows due to the mergers with 
more abundant NFW halos. This gives a conservative bound of the evolutionary effect of the halo density 
profiles. 
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60 I 1 
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50 - Zj = 1 .08 z, + 1 .2 




Figure 2: Initial redshift Zi of a iialo which has a probability 1 — 1/e when growing from initial 
mass M = 10^ Mq to 2M at zj. 



distribution of the final redshift p{zf,2M; Zi,M). The characteristic final redshift is defined 
as P{z > Zf,2M;zi,M) = 1 — 1/e, where P is the accumulative probability that the halo 



reaches 2M before Zf. Using the probability distribution Eq. (2.22) in Ref. |33] we find 
that an empirical fit Zi w l.OSzj + 1.2 nicely applies for the cosmological model adopted 
in this work, as shown in Fig. ^. As a result, the lower limit of redshift integration in Eq. 
( |3.5| ) is simply taken as l.OSz + 1.2. These two cases (infinite and minimum ages) bracket 
the possible scenarios of DS evolution history. 



Replacing Ltot in Eq. (3^) with L^^^ we can derive the enhancement factor A^^, for 



the halos with DS formation. The enhancement factors are shown in Fig. |^. It is shown 
that for the most optimistic case, i.e., without LW feedback and no evolution of the DS 
host halos, the enhancement factor due to DM dumpiness is found to be larger by 1 — 2 
orders of magnitude compared with the scenario without DS formation. However, for most 
other cases the additional enhancement effect due to DS formation is negligible. 

3.3 Gamma-ray flux and constraints on DM model parameters 

The 7-ray flux produced by DM annihilation observed today can be derived through inte- 



grating the emissivity over the evolution history of the Universe 
c nlp'iiav) ^ ,{l + z'f[l + A'^{z')]dN 



4tt 



2m^ 



dz' 



H(z' 



dE 



-exp [-t{z',E')] 



(3.6) 



where niy. is the mass of DM particle, {av) is the velocity weighted annihilation cross 
section, ~ 0.23 is the DM density parameter, pc = SHq/SttG is the critical density of 
the Universe at present, H{z) is the Hubble parameter, E' = E{l + z'), dN/dE' is the 7-ray 
generation multiplicity at redshift z' for one annihilation of a pair of DM particles, and 
t{z',E') is the optical depth of the 7-ray photons when propagating in the intergalactic 
space. 
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Figure 3: Enhancement factor of DM annihilation from dumpiness. Two concentration models, 
BOl (left) and power-law (right), are adopted. The red solid line is for the case without DS 
formation. Other lines are the enhancement factors for the halos with DS formed, for infinite 
(thick) and minimum (thin) age of the halos. From top to bottom the curves represent different 
LW efhciencies /lw = 0.0, 10"'', 10~^, 1.0 respectively. 



The 7-ray photons of both the primary component which are generated directly from 
the DM annihilation products and the secondary inverse-Compton (IC) component due to 
scattering of DM-induced electrons/positrons off the cosmic microwave background (CMB) 
photons are included. The photon spectrum of the primary component for specified an- 
nihilation mode is calculated using the simulation code PYTHIA |58|. For the secondary 



IC component, we assume the electrons and positrons will cool instantaneously after the 
production [59|. Then the equilibrium spectrum dne/dSg is simply the solution of the 



energy loss equation 



d 



~dfdK 



dNe 



where the energy loss rate dEjdt = 2.67 x IQ-^^l + z)'^{Ee/GeVf GeV s'^, dN^/dE^ 
is the production multiplicity per annihilation. Then the IC photon spectrum can be 
calculated through convolving the spectrum with the CMB photon spectrum and the 
Klein-Nishina differential cross section . 

The processes absorbing 7-ray photons include photo-ionization, photon-nuclei pair 
production, Compton scattering, photon-photon scattering and photon-photon pair pro- 
duction |]6T| , |6^ ]. For redshift lower than 6, we also consider the pair production when 
scattering off the cosmic infrared background, adopting the "baseline" model in Ref. [^]. 

We first consider the canonical neutralino-like DM model. The expected EGRB fluxes 
for = 100 GeV, (av) = 3 x 10~^^ cm^ s~^, bb final state are shown in the upper panels 
of Fig. ^. The left panel is for BOl concentration model, and the right one is for power-law 
concentration. The solid line denotes the ordinary model without DS formation. Several 
DS models are plotted to show the enhancement effect of DS formation: /lw = 0, T(Js = 00 
(long dashed); /lw = 10"'*, Tds = 00 (short dashed); /lw = with minimum r^s (dotted). 
We can see that for the case with the strongest DS enhancement the expected EGRB will 
exceed the Fermi data and should be excluded. On the other hand, for the models without 
DS formation or models with moderate DS enhancement, the data can not give very strong 
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Figure 4: Upper panels: comparison of tlie expected EGRB fluxes from DM annihilation -with 



the observational data by Fermi |18 . The particle parameters of DM are adopted as the canonical 



neutralino-like WIMP models -with = 100 GeV, (crw) = 3 x 10"^® cm~'^ s~^, and bb annihilation 
final state. Lo-wer panels: constraints on the — (av) parameter plane of the DM (regions 
belo-w the lines are allo-wed). The dots sho-w the random scan of the MSSM model parameters 
■with DarkSUSY [|4j , -with relic density of DM compatible -with (diamond) or lo-wer than (triangle) 
WMAP observational result. The left and right panels are for BOX and po-wer-la-w concentrations 
respectively. 



constraint. 

The exclusion limits on the rriy. — {av) plane, above which the parameter space is ex- 
cluded, are shown in the lower panels of Fig. ^. Here the constraints are derived by requir- 
ing that the EGRB not to exceed the 2a errorbars of the Fermi data, like the conservative 
method given in [ |65| . Actually since the spectral shape of Fermi data is very different from 
that expected from DM annihilation, we can assume a power-law background and expect to 
give a much stronger constraint on the DM contribution |^5|. For comparison we also give 
the supersymmetric model predicted parameters through a random scan of the parameter 
space of the Minimal Supersymmetric extension of the Standard Model (MSSM), using 
DarkSUSY package [^]. It is shown that for the model with the largest DS enhancement, 
a large part of the MSSM parameter space can be excluded. The exclusion limits even 
reach the range required by the relic density of DM if DM is produced thermally in the 
early Universe. 

What of most interest are the recent discoveries of the high energy electron/positron 
excesses by PAMELA H], ATIC Q, HESS p|, and Fermi 0. Together with the 
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Figure 5: Same as Fig. ^ but for leptonic DM model with /i^/i^ annihilation final state, which 
is proposed to be able to explain the recent observed electron/positron excesses. The mass of DM 
is = 1.7 TeV, and the annihilation cross section is (av) = 3.6 x 10~^^ cm~^ s~^. In the lower 
panels the contours are 3cr and Scr confidence regions of the fit to PAMELA/Fermi/HESS electron 
and positron data | |6^ . 



non-excess of antiproron-proton ratio 
responsible for the excesses (e.g., [38, 



, the leptonic DM model is favored if DM is 
7|, Inl, ^, 0, H, H, H). Alternatively, 



many astrophysical models have been proposed to interpret the signature (see [^] for a 
review). In this work we consider the DM annihilating to fJ-~^fi~ or t~^t~ final states, which 
are proposed to fit the PAMELA/Fermi/HESS data of positrons and electrons. 

The results for /^^^^ final state are presented in Fig. ^. In the upper panels, the 
illustrating EGRB fluxes for = 1.7 TeV and (av) = 3.6 x 10^^'^ cm^'^ s^^ are shown. 
The lower panels give the exclusion limits on the rriy. — {av) plane. The dashed contours in 
the lower panels are the 3a and 5a confidence regions of the fit to PAMELA/Fermi/HESS 
electron and positron data [js^ j. We can see that even for the case without DS enhancement, 
the favored parameter regions are excluded for power-law concentration model by the 
Fermi EGRB data. For BOl concentration model the exclusion limit is also very close 
to the favored region^ (note here the exclusion limits are very conservative). If the 
DS enhancement is non-negligible, the constraints will be even more stringent. For the 
DS scenario with strong enhancement (/lw = 0, t^s = oo) and moderate enhancement 



(/] 



LW 



10 
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= oo), the model explaining the e excesses can be well excluded. 



* Compared with the similar semi-analytical model "BulSub" of Ref. js^, the constraint is weaker here 
because the subhalos inside each halo are not considered. 
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Figure 6: Same as Fig. H but for leptonic DM model with t'^t annihilation final state. The 
mass of DM is — 3 TeV, and the annihilation cross section is (av) = 1.9 x 10^^^ cm^"^ s^^. 



The results for r+r^ final state are shown in Fig. ^. The DM model parameters chosen 
for illustration are = 3 TeV and (av) = 1.9 x 10^^^ cm^^ s^^. Similar conclusions as 
the fJ'^H~ case can be drawn. 



4. Conclusion and discussion 

In this work we study the constraints on the DM parameters using the Fermi measured 
isotropic 7-ray data, taking into account the formation of DSs accompanied with the Pop III 
stars. The DS formation is expected to result in a halo with enhanced density distribution 
of DM, which can give larger annihilation luminosity and produce higher 7-ray fluxes. Two 
kinds of DM particle models are discussed: the canonical neutralino-like particle, and the 
leptonic DM which might be responsible for the excesses. 

The formation rate of DSs is closely related to the SFR of Pop III stars. We employ 
an analytical way to calculate the SFR of Pop III stars. The LW feedback from the Pop III 
stars themselves and the galaxies is considered. However, there is large uncertainty of the 
self-shielding effect of the LW photons. A phenomenological LW efficiency parameter /lw 
is employed to parameterize different efficiency of the LW feedback. The results show that 
different /lw can give very different SFR of Pop III stars, and hence the enhancement effect 
of DM annihilation luminosity of the DS host halos is very different. For the case with 
little LW feedback (/lw ^ 10^"^) the enhancement due to DS formation can be remarkable. 
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Another large uncertainty is the evolution effect of the DM halos once DSs were formed. 
The halos will be subject to mergers during the evolution of the Universe. But the change 
of the density profile is not clear after the major merger of two halos with different initial 
profiles. Some of the numerical simulations suggest that the resulting density profile will be 
some intermediate profile after the major merger |54|, however, there are also simulations 
showing the final density profile will be essentially close to the cuspier one |^5|. Furthermore 
the accretion or minor merger may also affect the evolution of the halo profile. We adopt 
two extreme approaches to cover the evolution effects. One is that there is no evolution of 
the DS host halos after their formation. The other is that the DS host halos experience a 
fast evolution with minimum age, which means the cuspy density profile after DS formation 
will disappear soon due to one major merger. These two approaches will give very different 
enhancement effects of the DS host halos. We hope further numerical simulations will help 
to clarify this issue. 

There are other uncertainties such as mass function, parameters of first stars, density 
profile of DM halos after adiabatic contraction, fragmentation during the formation of the 
first stars. However, compared with the above two major uncertainties, these uncertainties 
are expected to be much smaller. For example, we tested that if the mass function is 
adopted as the PS form with {A, a, p) = (0.5, 1, 0) [^, the change of the enhancement 
factor is within a factor of 2. 

Given these large uncertainties, especially those from the SFR and halo evolution effect, 
our conclusion is also model dependent. We find that in the most optimistic case for 7-ray 
production through DM annihilation, i.e., weak LW feedback and no evolution, the DS 
host halos would enhance the 7-ray signals by 1 — 2 orders of magnitude. In this case, the 
constraints on the cross section of the neutralino-like DM can reach the thermal production 
range (~ 10~^^ cm^ s~^). The leptonic DM models proposed to explain the excesses can 
be well excluded, and the excess features would be most likely of an astrophysical origin 
(e.g., |77|). Conversely, if the excesses are of a DM annihilation origin, then the early 
Universe environment should be such that no DS can form. 

In other cases with stronger LW feedback and/or minimum age of the DS host halos, 
the DS enhancement is much smaller and the constraints are much weaker. However, for 
the power-law model, due to the large contribution of many low mass halos the constraints 
on the models to explain excesses still apply. 

We point out that our constraints are conservative to some extent in several aspects. 
First, the constraints are derived according to the 2a upper bounds of the error bars of 
the data. Since the observational spectrum is essentially a single power-law shape which is 
very different from the DM expected spectrum, one would expect the contribution of DM 
component to the EGRB would be even lower than the observational data |65|. Using that 
reduced background would give more stringent constraints on the DS formation scenarios 
and DM annihilation model. Furthermore, there are sub-halos and sub-sub-halos in each 
halo, which may give even stronger annihilation signals f?^. Considering this effect would 
pose even stringent constraints. Finally, only the adiabatic contraction processes accom- 
panied with the first star formation are included. There should be many other adiabatic 
contraction processes, e.g. during the formation of the second and third generation stars. 
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which should also give enhanced DM distribution. 
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